Optimal control of molecular electronics by femtosecond laser pulses 
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Combining the features of molecular wires and femtosecond laser pulses gives the unique oppor- 
tunity to optically switch electron currents in molecular devices with very high speed. Based on a 
weak-coupling approximation between wire and leads a quantum master equation for the population 
dynamics and the electric current through the molecular wire has been developed which allows for 
arbitrary time- dependent laser fields interacting with the wire. This formalism is combined with the 
theory of optimal control. For a tight-binding approximation of the wire we show how to compute 
the laser pulses to switch the current through the wire on and off. With this approach the desired 
pattern of the current in time can be chosen in an almost arbitrary fashion. 

PACS numbers: 73.63.Nm, 03.65.Yz, 33.80.-b 



The field of molecular electronics and especially of 
molecular wires has attracted much interest recently, ex- 
perimentally as well as theoretically P, • In scattering 
theories one usually uses the Landauer-Biittiker formula 
Q which was first developed for steady-state situations 
and later extended to deal with time-dependent fields 
[H HIi applying non-equilibrium Green's functions. In a 
different class of approaches one treats the coupling of the 
wire to the leads perturbatively [1, 0, H, H, EH EH ■ Some 
of these treatments utilize knowledge from the field of 
dissipative quantum dynamics and derive quantum mas- 
ter equations (QMEs) to describe the current through 
a molecular wire. But instead of coupling the system, 
i.e. the wire, as usual to a bosonic heat bath it is now 
coupled to a fermionic particle reservoir with which it 
can exchange particles. The obtained QMEs easily allow 
us to apply external optical fields which influence the 
dynamics in a direct manner by changing the wire part 
but also in an indirect way by influencing the wire-lead 
coupling. For monochromatic fields the Floquet theory 
has been employed @, El El El . We use a technique 
in the field of dissipative quantum dynamics in which a 
special parameterization of the so-called spectral density 
of the reservoir leads to a set of coupled equations for 
a primary and several auxiliary density matrices. This 
formalism allows for an arbitrary time dependence of the 
laser pulses [l5j . 

For cw-laser light it was shown [141 ] that, within a 
simple tight-binding model, the time-averaged current 
through a molecular wire could be suppressed choosing 
the right amplitude of the field. This is known as the phe- 
nomenon of coherent destruction of tunneling (CDT) and 
was first discovered in two-level systems [la ]. Recently 
we demonstrated [H EH that one can also suppress the 
current using short Gaussian laser pulses. The maximum 
amplitude of these pulses has to fulfill the same condi- 
tions as for a cw-laser pulse. Though the current can 
be suppressed the shape of the current pattern in time 
cannot be controlled easily. 

Lately an optimal control theory (OCT) for time- 



independent targets has been extended to control target 
states which are distributed in time [H H(| or which are 
time-dependent [21| . In this Letter we apply the OCT to 
molecular wires with a small number of sites. After se- 
lecting a time-dependent current target, one utilizes OCT 
to optimize the time-dependent control field. Under the 
influence of this field, indeed a current which is similar 
to the predefined current pattern is obtained. Although 
we concentrate on the topic of molecular wires here, the 
techniques described below can, in principle, be applied 
to coupled quantum dots as well. 

The time-dependent Hamiltonian of the investigated 
molecular junction is separated into the relevant system 
Hs(t), mimicking the wire, and reservoirs Hr modeling 
the leads 



H{t) = H s (t) + Hp(t) + H R + H t 



SR 



(1) 



with wire- lead coupling Hsr. Denoting the creation (an- 
nihilation) operator at site n by cj, (c n ) the tight-binding 
description of the electrons in the molecular wire reads 

Hs(t) = ^e„4c n - A^2(clc n+1 +4 l+1 c n ) . (2) 



The first term describes the on-site energies and the sec- 
ond term the nearest-neighbor hopping. For simplicity 
electron spin and interaction have been neglected since 
they will not qualitatively change the optimal control sce- 
nario. Denoting the dipole operator which will be de- 
tailed below by fi, the coupling between the wire and the 
laser field E(t) reads 



H F (t) = -fiE(t) 



(3) 



The environment of the wire consists of two electronic 
leads that are modeled by two independent reservoirs of 
uncorrelated electrons in thermal equilibrium. For each 
lead, the Hamiltonian Hr is given by Hr = Y^q^qbqbq 
with b q creating and b q annihilating an electron in the 
corresponding reservoir mode | q) with energy u> q . We set 



2 



h = 1. In further derivations we will only refer to the 
left lead but the formalism has to be applied as well to 
the right lead coupled to the last site N of the wire. The 
coupling of the left electronic lead with the first site of 
the wire is given by 

2 

Hsr = Yl K *** = E( + V g *bl Cl ) (4) 

x=l q 

with $! = J2 q V q b q , $ 2 = J2 q V*b\, K x = c\, K 2 = ei, 
and a wire-lead coupling strength V q for each reservoir 
mode. 

Starting with the time-convolutionless approach, a 
time-local QME based on a second-order perturbation 
theory in the molecule-lead coupling has been devel- 
oped for the reduced density matrix p(t) of the molecule 



oped toi 



dp® 

dt 



= -iC s (t)p{t) - iC FP (t) - D(t)p(t) , (5) 



D(t)p(t) = ^[K x k xx ,{t)p{t) - K xx ,{t)p{t)K x 

xx' 

-K x p(t)A xx , (t) + p{t)A xx , (t)K x ] (6) 



with auxiliary operators for the wire-lead coupling 



A xx >(t) 



dt'C xx/ (t-t')U s {t,t')K x , 



dt'C* xlx (t-t')U s (t,t')K x , 



(7) 
(8) 



Here we employed the definitions Us{t,t') — 

T + exp{- l / t !dr (C s (r)+C F (r))}, C s (t) = [H s (r),»], 

Cf{t) = [Hp (t),»] with the time-ordering operator 
T_|_ and the reservoir correlation functions C xx i(t). For 
later convenience we also introduce = \p, •]. Using 
the electron number operator Ni = c^ q c q of the left 
lead with the summation performed over the reservoir 
degrees of freedom yields 



Il(t) - e-tr {N tP (t)} = tr{I,(«)p(*)} (9) 



with the current operator 



2i(*)p(*) 



c\k 12 {t)p(t) - c\p(t)l l2 {t) 
c lP (t)A 21 (t) - ciAai(t)p(t)l (10) 



and the elementary charge e. This equation describes 
the current Ii(t) from the left lead into the molecule. 
A similar expression holds for I r (t) from the right lead 
into the molecule. As in a steady-state situation, the 
net transient current can be defined as I(t) — (Ii(t) — 
I r (t))/2. In addition, we define the net current operator 
l(t) = (l l (t)-l r (t))/2. 



Aim of this study is to determine laser pulses that re- 
sult in a predefined effect on the current through the 
molecular wire. This is achieved by extremizing a control 
functional. Using techniques developed previously for 
time-dependent targets [22], |23| we define as part of the 
functional the difference between a preselected current 
pattern P(t) and the current obtained from the QME 



J (E)= / dt{P(t)~ti[X(t) P (t)}} 2 



(11) 



To this we add a second part which shall ensure conver- 
gence SHi 



J(E) = ME) 



dt 



(E(t) - E(t) 



(12) 



with E(t) being the laser field of the previous iteration 
step. The penalty parameter A is a Lagrange multiplier. 
Introducing the time-dependent function s(t) we avoid a 
sudden switch-on and switch-off behavior of the control 
field at the beginning and the end of the propagation 
time. In the present study either a Gaussian or a squared 
sine function has been employed. The numerical results 
are rather independent of the actual choice. 

In order to calculate the extremum of the functional 
(TT2")) . the functional derivative has to vanish [22| : 



8J(E) = Q= 6Ja(E) 



A 6 



SE(t) 



5E(t) 2 6E(t) 



dt 



{E{t)-E(t)f 
s(t) 



(13) 

Using the time-evolution operator U(t,t',E) = 
T+exp j-i f* dr [C s {t) + C f {t) ~ i-D(r)]} together 
with p(t) = U(t 7 t 0l E)p(t a ), one gets 



E(t) 



E(t)- 
s(t) 



8® 

A 

(SO(t) 



A tT \SE(t) 



Pit) 



SD(t) 
5E{t) 



Pit) 



with the operators 

X (t) = J 1 drO(T)U(T,t,E), 
0(r) = 2{tr[J(r) /9 (r)]-P(r)}J(r) 



(14) 

(15) 
(16) 



The functional derivatives SD/5E and SO/SE are 
straightforward to evaluate but lengthy and will be given 
elsewhere. For the operator x(t) one can derive a QME 



dt 



-iC s (t) X (t) - iC F (t) X (t) + D(t) x (t) - 0(*jjl7) 



where the inhomogeneous term 0(t) reflects the time- 
distributed target state. The wire-lead coupling operator 
D(t) in this QME is given by 

D(t) X (t) = Y\ X (t)K x A xx ,(t) - K xX {t)A xx ,{t) 



-A xx > (t)x(t)K x + A xx > (t)K xX {t) (18) 



FIG. 1: Target current P(t) (dashed) and calculated current 
I(t) (solid) in the upper panel. The control field E(t) is shown 
in the lower panel. 



While Eq. §5§ is propagated forward in time, Eq. (JTTJ) is 
propagated backward in time. These two equations can 
now be solved in an iterative fashion. One starts with an 
equilibrated state on the wire p(0) at time to = which 
is unchanged during the iterative process and chooses an 
initial guess for the field E(t). The density matrix p{t) is 
propagated forward from initial time t = to to final time 
t = tf using E{t). In a second step, one obtains a new 
field E(t) from Eq. fl"4")) . in order to propagate the target 
operator x(t) backward from final time t — tf to initial 
time t = to- One stores the field E{t) as E(t) for the next 
step, in which we again calculate E(t) and propagate p(t) 
from to to tf with the initial condition p(to). This scheme 
is re pea ted iteratively until convergence. As discussed in 
Ref. [2l| some of the density matrices and fields can be 
stored during the iterative process to save computational 
effort. 

For the calculations shown below we consider a wire 
consisting of two sites without spin, equal site energies 
£1 = £2, intersite coupling A = 0.1 eV, and a dipole op- 
erator p — (c\ci—c\c2)/2. A bias Vt, = 0.4 eV symmetric 
with respect to the site energies e„ determines the Fermi 
energies Ef.i = £1 + V&/2 and ep, r = £2 — which 
describe the occupation of the reservoir modes. The cou- 
pling between the wire and the leads is almost in the 
wide-band limit and the maximum coupling strength is 

0. 1 A. For the first and third example below A was set 
to 0.003 and for the second calculation to 0.0006. 

As a first simple example of current control we con- 
sider a case in which the current is initially constant and 
the goal is to suppress the current following a Gaussian 
shape in time. In Fig. [T]the current target P(t) is shown 
together with the obtained current and the correspond- 
ing laser field E(t). Until the initial time to = the 
laser field is turned off and the system is equilibrated, 

1. e. in a steady state leading to a time-independent cur- 
rent. As can be seen the optimal control current matches 
the target very well. The optimal laser field is slightly 
asymmetric and its maximum occurs a few femtoseconds 
before the minimum of the current. This is in agree- 



FIG. 2: Same as in Fig. 1 for a different control target. 
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FIG. 3: A control target including a highly oscillating pattern. 
The control target (dashed) and the achieved current (solid) 
are displayed in the middle panel while in the top panel an 
average over five oscillatory periods is shown. The control 
field is given in the lowest part of the figure. 



ment with the earlier observation [171, [18| that a Gaus- 
sian shaped envelope of the laser field leads to a slightly 
asymmetric current pattern with a minimum delayed by 
a few femtoseconds. 

In the second example the complexity of the control 
task is increased. As shown in Fig. [2] the current target 
P(t) is a symmetric double step function and the goal is 
achieved by the optimal control algorithm rather accu- 
rately. Only at the last step there is some visible devia- 
tion between target and achieved current. The optimized 
laser field in the lower panel of Fig. O also shows step-like 
features but additionally oscillations and peaks. In con- 
trast to the previous control target, the rapidly changing 
target current pattern requires larger changes in the elec- 
tric field. To compensate the effect of these large peaks 
on the current at later times, additional oscillations of the 
field are needed to achieve constant values of the current 
at the plateaus. 

In many experimental setups the laser field has a high 
carrier frequency and a more slowly varying envelope 
function. To mimic such a scenario we created a target 
for the current which consists of a fast oscillating pattern 
with uj = 1 eV modulated by a slowly varying envelope 
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FIG. 4: The value of the control function J(E). Solid for the 
step function as shown in Fig.[2]and dashed for the symmetric 
Gaussian shape in Fig. [3] 

function, namely a current depletion of half-Gaussian 
form, a constant part, and again a half-Gaussian shape 
for the increasing current. This pattern together with 
the current achieved by the optimal control algorithm is 
shown in the middle panel of Fig. [3J Despite the rapidly 
varying target function the goal is accurately achieved. 
Similar to previous studies of CDT wc also calculated a 
mean current which is determined by averaging the cur- 
rent over a few cycles of the carrier frequency, here 5 
periods. This procedure can of course be performed for 
the target as well as for the shaped current and the re- 
sults are shown in the top panel of Fig. [3l As can be 
seen the average current is suppressed nicely in such a 
scenario. 

In previous studies [13, U3, CDT was utilized to 
suppress the current which is possible for certain values 



of the amplitude A of the laser field: The fraction A/lu 
has to be equal to a zero of the zeroth-order Bessel func- 
tion. In the present case the maximal value of the field 
is about A/uj = 0.8 which is certainly far off a zero of 
the zeroth-order Bessel function. Therefore the present 
current depletion cannot be explained by CDT. 

The convergence of the iterative control algorithm for 
the second and third example is displayed in Fig. [U It 
shows how the value of J decreases with increasing num- 
ber of iterations. For a simple target like a step function, 
J decreases very fast and then gets more or less constant. 
In contrast for a complex target with a fast oscillating 
function, it decreases not so fast in the beginning and 
later even increases slightly. To further decrease the con- 
trol function one could increase the Langrange parameter 
A for large iteration numbers. 

In conclusion, we have demonstrated that it is theo- 
retically possible to achieve a time-dependent pattern of 
the current through a molecular junction. The optimized 
current can be obtained using an optimal control field. 
Although we used a simplified model for a molecular wire 
the present investigation shows that it is worthwhile to 
study the combination of molecular wires and optimal 
control theory on the femtosecond time scale, theoreti- 
cally as well as experimentally. In experiment feedback 
coherent control for complex systems has been successful 
[26| so that an experimental implementation of the ideas 
described above might be feasible in the future. 

The authors would like to thank the Deutsche 
Forschungsgemeinschaft for financial support within SPP 
1243. 
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